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Abstract 

We provide a simple approach for the evaluation of inverse integral transforms that does not 
require any knowledge of complex analysis. The central idea behind the method is to reduce the 
inverse transform to the solution of an ordinary differential equation. We illustrate the utility of the 
approach by providing examples of the evaluation of transforms, without the use of tables. We also 
demonstrate how the method may be used to obtain a general representation of a function in the 
form of a series involving the Dirac-delta distribution and its derivatives, which has applications 
in quantum mechanics, semi-classical, and nuclear physics. 
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I. INTRODUCTION 



Integral transforms are usually presented to undergraduate physics majors well before 
they are exposed to them in any rigorous fashion in their mathematics courses. The primary 
motivation behind the integral transform is to re-write an ODE given in one representation 
(the original t-domain), to another representation (the s-domain, or image domain) where 
the solution is hopefully simpler; t and s are called conjugate variables. Once the solution 
in the s-domain is known, the inverse transform is applied to obtain the desired solution in 
the original t-representation. In this way, the integral transform is presented as yet another 
tool to add to the arsenal of techniques that may be applied to find the solution of linear 
ordinary differential equations. 

The usefulness of the integral transform method for the solution of an ODE (or a set 
of ODEs) is clearly predicated on the ability to perform the inverse transformation back 
to the original domain. However, if the inverse transform requires complex contour inte- 
gration methods (e.g., the inverse one-sided Laplace, two-sided Laplace, Mellin, and Weier- 
strass transforms), most undergraduate physics students will find themselves at an impasse; 
namely, students learning about methods for solving linear ODEs will likely have had no 
exposure to complex analysis, and therefore, be unable to perform the inverse transforma- 
tion without the use tables, or perhaps symbolic mathematics packages such as Maple© or 
Mathematica©. If the inverse transform cannot be obtained, the student is left with no 
choice but to attempt to solve the original problem using the set of techniques presented to 
them in their course(s) on linear ordinary differential equations. 

We wish to stress from the outset that the outcome of this tutorial does not change 
the situation described above; that is, if the problem is to find the solution to an ODE 
via an integral transform, and the inverse transform is not tabulated anywhere, then this 
article will be of little help. Rather, we have in mind a scenario in which a particular 
quantity of interest is explicitly formulated in terms of an inverse integral transform, and 
where the evaluation of the inverse transform ostensibly requires the use of complex contour 
integration methods. Our target demographic then, are undergraduate physics students with 
no exposure to complex analysis, but already familiar with the solution of linear ordinary 
differential equations. 

To this end, the rest of our paper is organized as follows. In the next section, we provide a 
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brief overview of general linear integral transforms, from which we focus our attention to the 
most common transforms arising in physics, where the inverse transform is a complex contour 
integral. We then present some simple rules which may be applied to reduce the complex 
contour integration problem to the solution of a linear ODE.- In Sec. Ill, we illustrate 
the method by providing several examples of inverse transforms involving complex contour 
integrals, which are either not easily found in standard references, or not readily obtained 
using symbolic mathematics packages. In all of our examples, the emphasis is on how to 
find the inverse transform, and not technical mathematical details, which we are confident 
will be covered in more formal mathematics courses. In Sec. IV, we demonstrate how our 
approach may also be applied to obtain a multipole expansion of a function, which would 
otherwise require a much more extensive mathematical analysis. Finally, in Sec. V we close 
the tutorial with some concluding remarks. 



II. INTEGRAL TRANSFORMS 

We consider a generic integral transform of the forrn^ 

f(s)=X[F(t)} = /" 2 K(s,t)F(t)dt , (1) 

where tGl, K(s,t) is the kernel, and f(s) the image of F(t) under the transform X. The 
inverse transform is given by 

F(t)=l- 1 [f(s)} = f 2 K- 1 (s 1 t)f(s) ds , (2) 

J si 



such that 



I 2 K- l (s,t)K(s,t')ds = 6(t-t'). (3) 

J si 

Equation ([3]) may be used to determine the inverse Kernel, K~ 1 (s,t) and, therefore whether 
the parameter, s, is real or complex. 

We now recall some properties of the integral transform. Linearity, X[F(t) + F(t)] = 
X[F{t)\ +X[G(t)], and multiplication by a constant X[AF{t)\ = AX[F(t)] follow from the 
integral definition of the transform given in Equation (JT|). Furthermore, it also follows 
immediately from Eq. (1) that 
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/*2 rti aiz 

(sK(s,t))F(t) dt; /«(*) = jf ^F(t)dt, 



(4) 



where the superscript in parenthesis denotes the order of the derivative with respect to the 
argument of the function. 

For the purposes of this paper, let us now focus on a class of Kernels (see also Sec. Ill C 
below), given by 

K(s,t) = Ce ast , a, C EC (5) 



which have the properties 



dK T ^ T ^ ldK 
atK , sK 



(6) 



ds a dt 

We thus observe that multiplication by s amounts to differentiating with respect to t and 
derivatives with respect to s amount to multiplying by t. It is a straightforward exercise to 
show that after n-fold applications of the multiplication and differentiation rules, we obtain 



A ( — A\ n 

s n f(s) = ^ + L 4 L X[F^(t)] , /W( a ) = an[t»F(t)] , 



where 



n 

A n = [^(-l) m (sa) n - m e sat F (m - 1) (t) 



m=l 



(7) 



(8) 



are surface terms resulting from n-applications of integration by parts. Kernels of the type 
in Eq. (5) also possess the so-called "shifting" property (let us take C — 1 for clarity), 



f * 



A 



a 



e ast [e At F(t)]dt = l[e At F(t)\ , 



At ■ 



and 



/oo poo 
e ast e~ As F(t)dt= / e a{t -^ )s F{t)dt = X 
-oo J —oo 



A 

F[t+- 

a 



(9) 
(10) 



In Table I, we present the various transforms that may be constructed by making par- 
ticular choices for the parameters discussed above. Table II provides a summary of some of 
the properties of integral transforms with Kernels of the form K(s,t) = Ce ast . 
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Transform Type 


K(s,t) 


h 


t 2 


K-\s,t) 


Si 


S2 




One-Sided Laplace 


e -st 





oo 


e st 

•2-Ki 


c — ioo 


c + ioo 




Two-Sided Laplace 


e -st 


— oo 


oo 


e st 

2-Ki 


c — ioo 


c + ioo 





Fourier 


e -ist 
V2^ 


— oo 


oo 


e iat 
V2^ 


— oo 


oo 






TABLE I: Integral transforms constructed after making particular choices for the parameters dis- 
cussed in the text. 



Property 


Image Function 


Inverse Integral Transform with K(s,t) = Ce ast 


I 


af(s) + bf(s) 


aF(t) + bG(t) 


II 


f(n) {s) 


a n t n F(t) 


III 


s n f(s) 


A n + ^W(f) 


IV 


e- As f(s) 


F(t + i) 


V 


f(s + i) 


e At F{t) 



TABLE II: A summary of some useful properties of integral transforms with Kernel K (s, t) = Ce ast 
A. The method 

We are now in a position to present the general method for reducing the evaluation of 
the inverse integral transform, to a linear ordinary differential equation. The method can 
be summarized as follows: 

1) Find the ODE or difference equation satisfied by the image, f(s), whose inverse 
transform, F(t), we desire. 

2) Term by term invert the equation obtained in step 1) using the properties listed in 
Table II above. 

3) Solve the resulting ODE/difference equation for F(t), thereby obtaining the desired 
inverse transform. 

It should now be clear why our method is not particularly useful for the study of 
ODEs. In the case of integral transforms applied to ODEs, one starts in the t-domain, 
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looking for the solution, F(t). The integral transform gives an algebraic equation for f(s), 
which if solved, is followed by an application of the inverse transform to obtain F(t) in the 
original domain. However, if we apply steps 1) and 2) above to f(s), the result will simply 
be a restatement of the original ODE in the t-domain, which is of no help, given that it was 
the ODE in the t-domain that motivated the transform to the s-domain in the first place. 

Instead, our method is designed to reformulate a problem which may not be easily ac- 
cessible to a student, e.g., finding the inverse transform via contour integration, to one for 
which the student is already well equipped to handle — the solution of a linear ODE. We 
are not suggesting that the solution of the resulting ODE will be simple (or even possible), 
but rather, the key point is that the student will be empowered to find the inverse transform 
using methods they are already comfortable, without having to rely solely on tables. In the 
next section we present several applications of the method, which will serve to both clarify 
its utility and applicability. 

III. APPLICATIONS 

A. Two-sided inverse Laplace Transform 

From Table I, the two-sided Laplace transform (LT) is given by 4 



It is assumed that the two-sided Laplace integral, Eq. (11), is convergent for Re s = c; 
the question of whether or not the integral may be convergent for other values of s may be 
answered by powerful theorems in Laplace transforms, which are not the focus of the present 
tutorial.- The inverse transform is given by 



The constant, c, must be greater than the largest real part of the zeros of the transform 
function. Therefore, for the two-sided inverse Laplace transform (ILT), the contour may be 
closed in either the left or right half plane, the choice being dependent on the function, F(t). 

The two-sided LT may in fact be used as the basis for all of the integral transforms 
presented in Table I. For example, if one puts s — > is in Eq. (11), the Fourier transform 




(11) 




c+ioo 



(12) 
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(FT) (aside from a constant) is immediately obtained.- The one-sided LT is likewise simply a 
special case of Eq. (11) provided we replace F(t) — > F(t)Q(t), where Q(t) is the unit Heavisde 
function. Moreover, the two-sided LT allows for the treatement of functions defined over 
t G R, in addition to loosening the restriction imposed by the FT, namely, that F(t) be 
absolutely integrable over t 6 R. It is also evident from Table I that the two-sided LT has 
very simple operational rules; that is, A n = 0, whereas for the one-sided LT, A n ^ and 
requires a careful application to the term- by-term inversion of the ODE satisfied by f(s) 
discussed above. 

1. A few known results 

To see the power of our approach in action, let us first demonstrate how to obtain some 
well known ILTs 4 by making use of the much simpler two-sided operational rules. 

Example 1. Consider the following image 

/(*) = ^ n >°- ( 13 ) 

The problem is to find the ILT of the image function, which by definition, requires the 
calculus of residues. However, if we apply our three-step procedure outlined above, we find 
that f(s) satisfies, 

S fU(s) + (n + l)f(s) = . (14) 
Using the properties listed in Table II, we see that the ODE for F(t) is given by 

-tF^{t) +nF{t) = , (15) 

which has the solution F(t) = Ct n . The constant C is determined by demanding that the 
two-sided LT of F(t) (a real single variable integral) is identical to the image function f(s). 
It is easy to see that we must have 

F{t) = B~ 

which is a well known result from the one-sided ILT tables. Note that our simple two 
sided operational rules naturally reproduce the one-sided results, and no complex contour 
integration has been used. 
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Q n+1 



t" 



T(n + 1 



-6(f) 



(16) 



Example 2. Next, consider the image function 

f(s) = s n , n>0, 



(17) 



whose ILT we wish to find using our ODE method. The image function satisfies the following 
first-order ODE, 

sfU(s)-nf(s) = 0. (18) 
Applying the same approach as in the previous example, we find that F(t) satisfies 

tF^(t) = -(n + l)F(t) . (19) 

Now, it is easy to show that F(t) = C/t n+1 (C is an integration constant) is a solution to 
Equation (19). However, 

/oo -St 
( 2 °) 

does not converge for any n > 0. Therefore, while we have obtained a valid solution to 
Eq. (19), it cannot be the two-sided ILT of f(s) since the forward transform, viz., Eq. (20) 
does not exist. If F(t) cannot be given by the solution of a function, what is the other 
possibility? As it happens, the only way to obtain a solution to Eq. (19), and have the 
forward transform exist, is to extend the class of functions over which the two-sided LT is 
defined to so-called tempered distributions.— Tempered distributions are beyond the scope of 
this tutorial, but fortunately, in the present case, the required distribution is the well-known 
Dirac-delta distribution. 

Using the fundamental properties for the derivative of the Dirac-delta distribution, we 
observe that 



e~ st 



S^ n \t)dt = s n . (21) 



We therefore conclude that 



F(t) = B- 1 [s n } = 5 {n) {t) . (22) 



We note that Eq. (22) is also a solution to Eq. (19), but in the sense of a distribution, 
meaning that both the left and right hand sides of Eq. (19) must be taken under an integral 
sign, so that F(t) = 5 (n) (t) acts on some arbitrary test function. 

The message to be taken from this analysis is that the existence of the forward transform 
must be ensured, over and above any formal solution to the ODE. This is a critical point, 
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as it serves to warn adopters of this method to be cautious in its application, since as we 
have just seen, the function F(t) = C/t n+1 , does not guarantee that the image, /(s), exists, 
in spite of being a formal solution to Equation (19). 

Example 3. Another useful, tabulated, image function is 

f(s) = exp(-ns) , n > . (23) 

The ODE satisfied by f(s) is given by 

/W( S )+n/(s) = 0, (24) 

and the ODE statisfied by F(t) is therefore 

- tF(t) + nF(t) = . (25) 

As in Example 2 above, the non-trivial solution to this equation can only be satisfied by the 
Dirac-distribution, viz., 

F(t) = # _1 [exp(-ns)] = 5(t - n) . (26) 
Example 4. Our final example for this subsection is to find the two-sided ILT of 

— k/s 

f(s) = - , k > . (27) 

s 

The ODE satisfied by f(s) is 

s 2 f w (s) + sf(s)-kf(s) = 0, (28) 

which implies that F(t) satisfies 

tF m (t) + F W (t) + kF(t) — . (29) 

It is an excellent exercise to show that the only acceptable solution to this second-order 
ODE is 

F(t) = J (2Vki)e(t) . (30) 

Once again, the simple two-sided operational rules have resulted in the correct one-sided 
inverse Laplace transform. 
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2. A few lesser known results 



Example 5. Consider 



f(s) = exp(as 3 



(31) 



where a is a real positive constant. This image is rarely found in tables of ILTs, but it is 
easily dealt with using our method. To wit, note that elementary differentiation yields 



fW( s )-3as 2 f(s) = 0. 

The term by term ILT of Eq. (32) gives 

F®(t) + ±-tF(t) = 0. 
3a 



(32) 



(33) 



The solution to the above second-order ODE is well-known to be given by a linear combi- 
nation of Airy functions, 



F(t) = CjAi 



-t 




-t 


+ C 2 Bi 


.(3a) 1/3 . 


_(3a) 1 /3_ 



(34) 



where C\ and C% are constants of integration. Since A n = for the two-sided LT, we require 
that 

-st tpU\ _ n (35) 



lim e- st F(t) = 

t— >±oo 



from which we find that C% = 0. The constant C\ is obtained by once again ensuring that 
the forward transform yields Eq. (31), with the final result being given by 



F(t) = B- L [exp(as d )} 



(3a) V3 



Ai 



(3a)V3 



(36) 



Example 6. We now wish to find the ILT for the following images 

Ms) = . 7 ! k f T r r \ = & esc as , (37) 

and 

f B (s) = — ff = -bcotbs , (38) 

where ks and T are constants, and b = nksT. Equations (37) and (38) have tremendous 
importance in the formulation of quantum mechanics, semi-classical^ and nuclear physics^ 



10 



at finite temperature, as will be discussed below. We begin with /f(s) and note that 
esc (a; — 7r) = — esc x , which allows us to write 

/f (a - |) =-/*■(«)■ (39) 
We next multiply both sides of Eq. (39) by s — ? to obtain 



after which a term by term ILT yields 



(ex + l)F«(t) = ^(t) . (41) 



The solution to Eq. (41) is formally given by 



7T 



dt 



It is straightforward to show that 

Ff ® = e -Wb + I = e -t/k B T + I ' ( 43 ) 

where C is a constant of integration. We may fix C by appealing to final value theorem for 
the LT, namely, 4 

lim F(t) = \imsf(s) , (44) 



which in the present case gives 



e c = lim — i— = 1 . (45) 
s^o sm(os) 



Finally, it follows that 



F F (t) = B- 1 



1 (46) 



e -t/k B T + ^ • 



sin(7rfcsTs) 

We have only seen a similar form of Eq. (46) explicitly tabulated in one reference, where it 
is simply stated without derivation.- 
We now move onto 

/*(*) = " I nk f T ■ (47) 
tan it Kb I s 

Making use of the fact that cot (a; — n) = cot a;, we readily obtain 

Sb { s ~l)= Sb{s) , (48) 
11 



which after multiplication on both sides by s — ? yields 

fB(s-^) = ( S - 7 UMs). (19) 



The term by term ILT of Eq. (49) gives 

(e^ b -l)F^(t) = -^F B (t). (50) 
Analogous to the treatment of Fp(t) above, the solution to this ODE is given by 

^HCexp^/^).-^. (51, 
Another application of the final value theorem shows that (7=1 and therefore^ 



F B (t) = B- 1 



tan(7r ksTs) 



-t/k B T _ I ■ 



(52) 



3. Physical Application: Ideal harmonically trapped quantum gases at finite temperature 

An immediate application of the inverse transforms derived above are harmonically 
trapped, ideal quantum gases at finite temperature. A self-contained discussion of the 
physics of these systems is well beyond the scope of this article, so here, we wish to simply 
provide the basic mathematical tools, and introduce how the method we have proposed can 
be used to introduce a complementary formulation of quantum systems, that is not widely 
known to many physicists. The interested reader is encouraged to examine Refs. 8 . 9 ] and 
the references therein for more details. 

The central theoretical tool from which all thermodynamic properties of the ideal gas can 
be derived is the Bloch density matrix (BDM). At zero temperature, the BDM is given by 

C (r,r';s) = ^ V*(r)Vi(r') expose,) , (53) 

all i 

where the ^'s and the one-particle energies 6j are the solutions to the time-independent 
Schrodinger equation, and the quantity, s, is a complex variable. The BDM appears to 
require the explicit single-particle wave functions for its evaluation, but in fact, this is not 
always the case. The BDM satisfies the Bloch equation 

H r C (r, r ; s) = — , (54) 
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where H r is the Hamiltonian and Co(r, r'; s) is subject to the initial condition 

C (r,r';s = 0)=5(r-r') . 



(55) 



Therefore, an explicit form for Co is, at least in principle, possible once the external potential 
is defined, and used in the Hamiltonian, H r , to solve Equation (54). 

In the present application, we will take the external confinement to be a three-dimensional 



isotropic harmonic oscillator, V(r) = ^mu^r 2 , with r = x 2 + y 2 + z 2 , and Uq the trapping 
frequency. Under such a confining potential, the exact zero temperature BDM is found to 
be given by 8 



C (r,r'; S ) 



2-7T sinh(s) 



3/2 



exp 



r + r' 



tanh(s/2) 



coth(s/2) 



• (56) 



The zero temperature BDM is independent of the quantum statistics of the particles, mean- 
ing that we have a unified treatment for either Fermi or Bose gases (including Bose-Einstein 
Condensation). Using Eq. (56), we may obtain the finite temperature first-order density 
matrix (FDM) via a two-sided ILT,-&2 



p(r,r';T) = B- 1 |c T (r, r'; s) 



(57) 



where 



C T (r,r';s) = C (r,r';s)- 



-Cn(r, r; s 



sm(irsk B T) 



(fermions) 
(bosons) , 



(58) 



tan(7rs/csT) 

and now ks and T take on the physical significance of the Boltzmann constant and tem- 
perature, respectively. Note that in Eq. (57), we have introduced the subscript, /i, in the 
two-sided ILT to emphasize that the conjugate variable to s in this problem is /x, the chemical 
potential. The factor, g in Eq. (57) depends on the spins of the particles through g = 2si + l, 
where s\ is the spin of the particle. For spin-1/2 fermions, g = 2, and for spinless bosons, 
g = 1. The finite temperature BDM, CV, encodes the quantum statistics by appropriately 
weighting the zero temperature Bloch density matrix. It should now be evident why we 
must utilize the two-sided ILT, rather than the one-sided ILT; the chemical potential, /i, 
must be allowed to take on both positive and negative values at finite temperature, whereas 
the one-sided ILT only allows for \i > 0. 
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Equation (57) is the fundamental quantity from which all of the thermodynamic proper- 
ties of the trapped gases (fermions or bosons) may be investigated. For example, the finite 
temperature spatial density, p(r;T) is immediately obtained by putting r = r' in Eq. (57), 
while the the finite temperature kinetic energy density can be found from 

3 r 



r r 



;T> = E 



i=l 



o 2 



dxidx' A 



7 p{r, r; T) 



(59) 



Following the calculations in Ref. [9j, the ILTs we have presented earlier may be used to 
explicitly evaluate p(r, r';T). 

In the case where the potential, V(r), is such that an exact, closed form expression for 
Cq(t, r'; s) is not possible, an ^-expansion may be performed, viz., 



C (r,r';s) 



-I e- sV ^e 



2h 2 s 



h 2 s 2 

12m 



VV 



2 {VV)2 . 



(60) 



where for economy of space, we have only shown terms up to relative order ti 2 , A = 
(27th 2 s/m) 1 ' 2 and h is Planck's constant divided by 2n. Equation (60) is called the semi- 
classical expansion of the BDM, and the higher order ^,-terms, containing gradients of the 
potential, represent quantum corrections to the classical result. The semi-classical expansion 
may be used as the foundation for the study of quantum effects in a variety of systems by 
systematically introducing higher-order quantum corrections and examing their influence on 
the physical properties of the system. A detailed discussion of the applications of Eq. (60) 
in semi-classical physics may be found in Reference js]. 



B. Inverse Fourier transform 



The inverse Fourier transform (IFT) is mathematically much easier than the ILT, since 
the IFT is an integral over the real line. Nevertheless, it is worthwhile illustrating our 
method as it applies to the IFT, particularly in cases where standard approaches to finding 
the inverse transform require some "tricks" for its evaluation. 



Example 7. Let us then consider the following image 



./(.s) = A /^. (61) 



14 



Simple differentiation shows that 



r( l), , / /7T ^ \ / s \ s 



m . )= y-eTT (62) 

and thus f(s) satisfies the following first order differential equation 

f W (s) + ^f(s) = 0. (63) 

Applying the properties in Table II to Eq. (63) above leads to the following ODE for F(t) 

F {1) (t) + 2atF(t) = . (64) 

Obtaining the solution to Eq. (64) is elementary, and is given by 

F(t) = Ce- at2 . (65) 

The constant C is adjusted such that the Fourier transform of F(t) is f(s). It then follows 
that 



F(t) = F- 1 



IT ^sf 

— e 4ct 



e- at . (66) 



Obtaining Eq. (66) using standard techniques is more involved.— 
C. Mellin transform 

The Mellin transform (MT) may also be treated by the ODE method presented here. 
The Mellin transform is defined by^ 

POO 

f(s) = M[F(t)}= / t s - l F(t)dt, (67) 
Jo 

and its inverse is given by 

-r i-c+ioo 

F(t) = M~ l [f(s)} = — / t-/00 da . (68) 

Am j c _j 00 

Once again, we observe that the two-sided LT may be viewed as the fundamental tool for 
our operational calculus, since the MT may also be obtained from to the two-sided LT by 
noting that^ 

M[F(t)] = B[F(e- t )} . (69) 
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The only thing left to do is to determine the corresponding "rules" as we did for transforms 
with exponential kernels. It is straightforward to confirm that following rules are obeyed for 
the Mellin transform: 



M[F(at)\ = a- s f(s) 

M[t a F{t)} = f(s + a) 

M[FW(t)] = -(s-l)f(s-l) 

M[t n F^(t)] = (_i)»E^f)/( a ) . (70) 

Example 8. To illustrate the method applied to the MT, let us use it to evaluate the 
complex contour integral, 

^L r * r(s)ds - (7l) 

which we recognize as the inverse Mellin transform of f(s) = T(s). According to the rules 
given above, we note that 

f(s + l) = sf(s) . (72) 
Applying the inverse rules to Eq. (72) yields 

F W (t) = -F(t) , (73) 

which is immediately seen to to have the solution 

F(t) = Ce~ l . (74) 

The constant, C, is found by ensuring that the forward transform, f(s) = T(s) (which is just 
a real integral), is obtained, which gives C — 1. Therefore, we have evaluated the complex 
contour integral, 



- / t-V(s) ds , (75) 

J c—ioo 

using only elementary methods. 
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IV. MULTIPOLE EXPANSIONS OF FUNCTIONS 

In this section, we wish to illustrate how to obtain the "multipole" expansion of a func- 
tion, Fit), in the form of a series involving the Dirac-delta distribution and its derivatives. 
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Such expansions have been found to be useful in semi-classical physics, where they may be 
used to introduce quantum corrections,— r — and in nuclear physics, for the calculation of 
approximate wave functions for realistic nuclear potentials.— 

The method presented here is general, and uses a very simple approach based on the two- 
sided inverse Laplace transform. An immediate benefit of our formulation is that it permits 
an easy determination of all of the moments of a function, F(t) (see Eq. (79) below). In 
addition, our approach allows us to obtain results that have previously required intense 
mathematical analysis and dedicated papers for their derivation. 

In the most generic of cases, we begin with the forward two-sided LT of F(t), viz., 



oo 

St 



f(s) = / e- st F(t)dt . (76) 

J — oo 

Assuming that f(s) is analytic, we use the series expansion, e~ st = ^2^ =0 {—^) n s n t n /n\, in 
Eq. (76) to obtain 

/OO = E f^T r tnF (t)dt) s». (77) 

n=0 V n - J -°° ' 

At this point, we are not worrying about the fact that we have taken the summation outside 
of the integral, since our end result will be an expansion that is defined only in the sense 
of a distribution, and therefore will only be valid under the integral sign. We now term by 
term Laplace invert Eq. (77), which yields 

00 / -i\n 

F(t) = Y, [ -^ L U {n) (t) , (78) 

n=0 

where 

/oo 
q n F(q) dq , (79) 
■oo 

is the n-th moment of the function, F(t). The n = term is sometimes referred to as the 
the "monopole" moment, and the n— \ the "dipole" moment, and so on. Equation (78) is 
the general form of the multipole moment expansion, and is to be viewed in the distribution 
sense, having meaning only under the integral sign. 

Example 9. In most practical applications, the moments of the function, Eq. (79), 
need not be explicitly calculated. To see this, let us first consider F(t) = Ai(t). We can 
deduce from Eq. (36) that, 

f(s) = exp(^-\ , (80) 
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is the two-sided LT of F(t). We then formaly write f(s) as a series expansion, 

v 7 n=0 



and then term by term ILT both sides of Eq. (81) (with the aid of Example 2) to obtain 
Ai(t) = B- 1 



00 ( l\n e 3n 



3 n n\ 

n=0 



■. y tirs^M . ( 82) 

^ 3"n! v ; 

n=0 

A comparison of Eq. (79) with Eq. (82) may be used to determine the moments of the Airy 
function. Equation (82) has been derived earlier by several authors,— r — but by using a 
much more involved approach than the one presented here. As commented above, it is im- 
portant to remember that Eq. (82) is to be interpreted strictly in the sense of a distribution. 

Example 10. Next we find an analogous expansion for the Bessel function, 

F(t) = J u (ojt)Q(t), where Re [v] > —1. The well known Laplace transform of this 

function is^ 

/(*) = B[J v (vt)e(t)] = (V " 2 t^Lj ■ (83) 

As in the previous example, we consider a formal series expansion of f(s), 

w"Vw 2 + s 2 ^ n!r (^f^) ' 1 ' 

followed by a term by term two-sided ILT of both sides of Equation (84). The result of this 



inversion is 

Uut)Q{t) = y 2 J° [ . (85) 

n=0 \ 2 / 



To our knowledge, Eq. (85) is a new result, and its use can immediately be seen from an 
examination of integrals for the form 



I{v) = / (p(r)J v (r) dr 



oo 2»-i r (5^)^(0) 

U ] nir(=^) 



(86) 



Again, a comparison of Eq. (79) with Eq. (85) allows us to obtain all of the moments of J u , 
viz., 
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which by direct application of Eq. (79) would have been more difficult. 



Example 11. As our final example, we consider the multipole expansion of the Fermi 
function,— 

f M - [e (,-H>/. + 1] ' I 88 ' 

In nuclear physics applications, R would denote the nuclear radius, and a the "diffuseness" 
of its surface. It is straightforward to deduce from Eq. (46), and property IV of Table II, 
that F(t) has the image 

f(s) = -^—,e~«* . (89) 

We next expand f(s) as a power series about the diffuseness parameter, a, and obtain the 
expansion 

— Rs 00 

" " - - ' - A 2k+1 a 2k+2 e- Rs s 2k+1 , (90) 



sin(7ras) 
where 

A - ~ ^^ IH I f91l 
Mk+l ~ (2TT2)! ' (91) 

and the B 2 k+2 are the known Bernoulli numbers. The term by term two-sided ILT of Eq. (90) 
gives 

^ oo 

[e( t-*)/ a + x] = ®( R - - E M k+l a^H^\t - R) . (92) 

Equation (92), which we have derived in only a few lines, has previously been the focus of 
an entire paper by L. Lukyanov,— thereby emphasizing just how useful the approach we 
present here can be in a wide array of problems. In fact, what is remarkable is just how 
trivial Eq. (92) was to obtain whe, oo.pated to the a„a lyS is p re8 e„ted in Reference Q. 



V. CLOSING REMARKS 



We have presented an alternative approach for the evaluation of common inverse trans- 
forms occurring in physics without having to introduce complex contour integration. Our 
operational calculus has been based on the two-sided Laplace transform, which has allowed 
us to examine a much wider class of functions than would be possible using other trans- 
forms. We have also highlighted the power of the two-sided inverse Laplace transform by 



19 



demonstrating how it can often make trivial work of problems (e.g., Examples 9 and 11), 
which otherwise require a much more extensive analysis. 

Some may find the presentation in this tutorial lacking in mathematical rigor. In partic- 
ular, we have purposefully avoided technical discussions of strips of convergence, for which 
at best, the integral transform is valid. For the examples we have provided, where the strip 
of convergence is not indicated, the latter has to be investigated, given the original F(t), 
for each individual case separately. Indeed, a useful exercise to give to students would be to 
determine the strips of convergence for each of the examples provided in this paper. 

Nevertheless, it is important to keep in mind that in undergraduate physics culture, 
advanced mathematical tools (e.g., integral transforms) are often needed right away — stu- 
dents just need to know how to use the tools, with the details of all the theorems and proofs 
being covered in more formal mathematics courses. Therefore, we hope the techniques pre- 
sented in this tutorial will also find their way into the physics classroom, where they may be 
used to enable students to evaluate integral transforms with confidence, and without being 
slaves to tables. We also hope that our presentation will inspire students to further examine 
the formulations of quantum mechanics, semi-classical physics, 6 and density-functional the- 
ory,-*^ all of which may be based on the two-sided Laplace transform. Finally, we believe 
that our general results for the multipole expansion of functions will be of interest to prac- 
ticing physicists, particularly those involved in the foundations of quantum, semi-classical 
and nuclear physics research. 
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After the completion of this work, we became aware of a long forgotten paper (Ref. [2j) in 
which a similar idea is proposed. However, Ref. [3] focuses soley on the one-sided ILT, which 
amounts to a special case of the present paper. Indeed, our presentation is applicable to a 
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wider array of functions, has simpler operational rules, and thereby allows us to present new 
results (e.g., the multipole expansions presented in Examples 9-11, which in our opinion, are 
publishable results in their own right) which would not be possible using the method presented 
in Reference Qj. Regardless, given that the approach presented in this tutorial is not discussed 
in any undergraduate textbooks we are aware of, it is certainly worthy of presenting to a broader 
audience. 
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